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ABSTRACT 



Context. We analyse the NGC 2024 Hn region and molecular cloud interface using [ 12 Cn] and [ 13 Cn] observations. 

Aims. We attempt to gain insight into the physical structure of the interface layer between the molecular cloud and the Hn region. 

Methods. Observations of [ 12 Cn] and [ 13 Cn] emission at 158 fim with high spatial and spectral resolution allow us to study the detailed 

structure of the ionization front and estimate the column densities and temperatures of the ionized carbon layer in the PDR. 

Results. The [ 12 Cn] emission closely follows the distribution of the 8 /im continuum. Across most of the source, the spectral lines 

have two velocity peaks similar to lines of rare CO isotopes. The [ 13 Cn] emission is detected near the edge-on ionization front. It 

has only a single velocity component, which implies that the [ 12 Cn] line shape is caused by self-absorption. An anomalous hyperfine 

line-intensity ratio observed in [ 13 Cn] cannot yet be explained. 

Conclusions. Our analysis of the two isotopes results in a total column density of N(H) ~ 1.6xl0 23 cm -2 in the gas emitting the [Cn] 

line. A large fraction of this gas has to be at a temperature of several hundred K. The self- absorption is caused by a cooler (T<100 K) 

foreground component containing a column density of N(H) x 10 22 cm -2 . 

Key words. ISM: atoms - ISM: clouds - ISM: individual objects: NGC 2024 - photon-dominated region (PDR) 



1. Introduction 

NGC 202 4 is a well- studied star forming region at a distance 
of 415 pc ( Ant hony-Twarogl 1 982) in the O rion B complex. Th e 
source consists of a dense (n « 10 6 cm -3 , Schul z~et~aD (Il99ll) ). 
narrow (~1'), north- s outh extended molecular cloud with an 
embedded Hn regi on (iKriigelet al.lll982h ICrutcheret all 1 19861: 
iBarnes et"aDll989l) . which extends several arcminutes beyond 

■ the molecular ridge in the east-west direction. The molecular 
material in front of the Hn reg ion is seen as a p rominent dust 

! lane in the optical. B oth OH (IBarnes et al.1 119891) and H 2 CO 
(iCrutcher et al.l fl986) absorption line measurements relative to 

| the radio continuum indicate that this material is at a radial ve- 
locity of vlsr ~ 9 km/s. The bulk of the molecular gas behind 
the Hn region is found at vlsr ~ 1 1 km/s from observations 
of e.g. optically thin CO (Mauersbergeret al. 1992; Grafetal. 
11993b Emprechti nger et all 120091: iBuckle et al,ll2010b or HCO + 
(JRicheret al.ll 19891) lines. In the south of the Hn region, a sharp 
ionization front delineates the boundary between the ionized 
and th e molecular gas (e.g. MSX data in Watanabe & Mitchell 
(2008)). This ionization front is seen edge-on in the south of the 
Hn region, but probably also extends north in a face-on geom- 
etry along both the foreground and the background parts of the 
cloud. 

The velocity dispersion in each of the two main parts of 
the cloud is small Avlsr ~ 2 km/s, which leads to a con- 
spicuous double-peaked line shape in many molecular emis- 



sion lines. Detailed modeling of these lines dGraf et al.l [l 993: 
Emprechtinger et al. 2009) results in H 2 column densities of 



a few times 10 22 cm" z for the foreground component and 1- 
2xl0 23 cm -2 for the background component. Gas temperatures 
are found to be around 70 K in the background and around 30 
K in the foreground, with a likely increase near the Hn region 
interface. 

In this Letter, we present [Cn] maps of an area near the ion- 
ization front. Our analysis concentrates on the physical proper- 
ties implied by the strong [ 13 Cn] emission detected. 



2. Observations 

We observed the 2 P3/2 — > 2 P\/2 fine structure tra nsition of ion- 
ized c arbon (C + ) at 1900.5369 GHz (157.7 /mi) dCooksv et al.l 
119861) toward the NGC 2024 (Orion B) molecular cloud. The 
measurements were made on Nov. 8, 2011 using the German 
REceiver for A stronomy at Terahertz frequencies (GREAT)Q. 
(JHeyminck et al .120121) on the Stratospheric Observatory for Far- 
Infrared Astronomy (SOFIA) (lYoung et al.ll2012h . 

The receiver noise temperature during the observations was 
approximately 4000 K (SSB). With a zenith opacity of 0.03- 
0.07, the system temperature was lower than 4500 K (SSB). The 



1 GREAT is a development by the MPI fiir Radioastronomie 
and KOSMA/Universitat zu Koln, in cooperation with the MPI fiir 
Sonnensystemforschung and the DLR Institut fiir Planetenforschung 
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beam size was 16". The velocity resolution before smoothing 
was 33 m/s. 

We mapped a 192" x 144" area in total power on-the-fly 
mode. Both the sampling distance along the scan direction and 
the scan spacing were 6" . The map was centered on (a, S) = 
(05 h 41 m 45 s .20, -01°55'45'/0) (J2000) using an off-pos it ion at 
(05 h 41 m 36 s .30, -02°05'0070). According to UaffeetaU (119941) 
and confirmed by a comparison measurement to a second, far- 
away off-position, no [Cn] emission could be detected at this 
position. 

The measurements were calibrated using the procedure out- 
lined in Guan et al. (2012). Instead of the main beam efficiency 
of 0.51 dHeyminck et al .112012b . we used a value of 0.58 to ac- 
count for the non-negligible coupling of the first side lobes to 
an extended source, as derived from a model of the diffraction 
pattern. A linear baseline was subtracted from each spectrum. 



3. Results 

Figure \T\ shows the measured [Cn] intensity distribution inte- 
grated over the velocity interval 5-15 km/s. T he emission agrees 
remarkably well with the 8 /mi MSX data dPriceetalJl2001t 
IWatanabe & Mitchellll2008h . which traces the UV-heated dust in 
the interface between the Hn region and the molecular cloud. 
Only the embedded MSX source just north of the map center is 
not very prominent in the [Cn] data. The emission peaks along 
the southern rim of the Hn region, where the interface is seen 
edge-on in a sharp ionization front. 

The intensity distribution does not reflect the distribution of 
the molecular material located in a roughly north-south elon- 
gated molecular ridge - outlined by the FIR peaks in Fig. Q] 

Averaged over a 60" diameter, we obtain a peak surface 
brightness of 4.4xl0~ 3 ergs cm " 2 s" 1 sr" 1 , in reasonable agree- 
ment with KAO and ISO data JJaffe et alJll99l iGiannini et all 
2000). 
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3.1. [ 13 Cn] distribution 

With a velocity span of ±100 km/s, the spectra also cover 
the three hyperfine structure components of the [ 13 Cn] transi- 
tion (vLSR-offsets -65.2, 11.2, 63.2 km/s dCooksvet al.lll986l) ) 
(Fig. [3]). If we co-add the velocity ranges [-60,-50], [19,25], 
and [70,80] km/s (taking into account the systemic velocity), we 
find emission with a similar spatial distribution as the main iso- 
tope (Fig.[T]). 

The [ 13 Cn] emission is only prominent in a crescent- shaped 
area following the edge-on ionization front, where column den- 
sities are highest. The peak emission coincides with the [ 12 Cn] 
emission maximum, but the [ 13 Cn] intensity falls off more 
rapidly, both to the Hn region side (north) and to the molecu- 
lar cloud side (south). The FWHM of the [ 13 Cn] emission zone 
is «!', corresponding to a linear size of ^0.1 pc. 



3.2. [ 12 Cn] channel maps 

The [Cn] emission is separated into two distinct velocity com- 
ponents, which show up as two peaks in most spectra (Fig. [3]). 
Fig. [21 presents channel maps of 1 km/s wide velocity bins rang- 
ing from 6 km/s to 14 km/s. The maps clearly show the two main 
velocity peaks, which are predominant throughout most of the 
source. The low velocity peak dominates the 8-9 km/s panels, 
the high velocity peak being clearly visible at 11 and 12 km/s. 
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Fig. 1. Top: Integrated (5-15 km/s) [Cn] intensity map of NGC 
2024 (color code). Comparison data: MSX Ban d A (8.28 //m) 
(con tours), 6 cm continuum peaks (triangles) (iCrutcher et al.l 
1986J), 1.3 mm c ontinuum emission clumps FIR2 to FIR6 
(squares) (Mezger et al. 1988), and embed ded infrared point 
sources IRS2 (lGrasdalenlll974 and IRS3 (iBarnes et alJll989b 
(stars). Bottom: Overlay of the [ 13 Cn] (color coded) and the 
[ 12 Cn] integrated intensities (contours: 100 K km/s to 625 K 
km/s in steps of 75 K km/s). The [ 13 Cn] map has been smoothed 
to a resolution of 25". The dashed box outlines the area of strong 
emission that was used in the analysis of section |4] 



The local minimum in the 10 km/s velocity bin around FIR5 co- 
incides with the peak of the integrated intensity map (Fig.[T]). 

There appears to be an anticorrelation between the low ve- 
locity peak and the high velocity peak. The relatively high in- 
tensities found at 9 km/s west of FIR2 to FIR4 and east of FIR5 
correspond to areas of weak emission at 12-13 km/s. 



4. Discussion 

For the data analysis, the interpretation of the two velocity peaks 
of the spectra is crucial. Although the [Cn] peaks - at velocities 
of ~8.4 km/s and ~12 km/s - are slightly offset from the two 
emission components in the optica lly thin lines of ve ry rare CO 
isotopes (9.2 km/s and 11.2 km/s. iGraf et al.l (1 19931) ). it seems 
reasonable to assume that the [Cn] line may also be caused by 
two emission components and not by self- absorption. The [ 13 Cn] 
measurement allows us to clearly distinguish between these two 
possibilities. 
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4.1. [ 13 Cn] spectrum 

To obtain a high signal-to-noise spectrum to analyze the isotopic 
line, we averaged the spectra within a rectangular area encom- 
passing the peak of the [ 13 Cn] emission (marked by a dashed box 
in Fig. []]). The resulting spectrum is shown in Fig. [3] In Fig. |U 
we compare the line profile of the [ 12 Cn] line to the three hyper- 
fine satellites of [ 13 Cn]. 

It is immediately obvious that the isotopic line only shows 
a single peak. Thus, either a) only one of the two velocity 
components seen in [ 12 Cn] has a sufficiently high column den- 
sity to produce measurable [ 13 Cn] emission, or b) the shape 
of the [ 12 Cn] is produced by a warm broad background com- 
ponent shadowed by a narrow cooler foreground component. 
Explanation a) would imply that there is an error in the rest 
freque ncies of ^10 MHz, clearly beyond the confidence lim- 
its of ICooksy et al.l dl986l) and in contrast to the good agree- 
ment found in recent H IFI measurements toward the Orion bar 
(lOssenkopf et al.ll201lh . Thus, if we accept the rest frequencies 
to their published accuracy, we are forced to the conclusion 
that the [ 12 Cn] line profile is caused by self-absorbed strong 
background emission centered at v L sr ~ 10 km/s. 

Self- absorption in the [ 12 Cn] line offers a quite natural ex- 
planation of the anticorrelation observed in the spatial distribu- 
tion of the two velocity peaks: A velocity shift of the absorbing 
component with respect to the background component shifts the 
observed line intensity from one velocity peak to the other. 

That the highest contrast between the 10 km/s dip and the 8 
km/s and 12 km/s pe aks is found in an area of kno wn strong fore- 
ground absorption (|Cmtcheretai][l986l iBarneseFal. 1989) also 
stren gthens this interpret ation. Radio recombination lines of car- 
bon (iKriigelet al .11 19821 and reference therein) are also mostly 
near 10 km/s. 

The discrepancy between the CO velocities and the self- 
absorbed [Cn] emission indicates that the simple source model 
with only two velocity components may underestimate the kine- 
matic complexity of the source. 



4.2. Hyperfine line ratios 

Figure |4] includes the result of a five-component Gaussian fit to 
the spectrum. The first two components model the main isotope 
(fit not shown in Fig. [4]) as an emissive and an absorptive com- 
ponent. The remaining three components fit the three hyperfine 
lines, where the line width is kept fixed and equal to the width 
of the broad background component. Table [T] lists the fit results 
and compares the parameters of the [ 13 Cn] hyperfine satellites to 
the literature values. 

For the two stronger hyperfine components, the velocities are 
within 0.5 km/s (3 MHz) of the predicted values. The signal-to- 
noise ratio of the F = 1 — » 1 component is too low to make a 
firm statement. The predicted line strength ratio, however, is not 
reproduced by our measurement. The (F = 2 — > 1)1 (F = 1 — > 0) 
intensity ratio should be 1.25, bu t is found to be 1.9. A similar 
discrepancy was also found by lOssenkopf et al.l (120111) . 

Obviously it cannot be ruled out that an as yet unknown 
emission line may coincide with the F = 2 —> 1 line, lead- 
ing to excess emission at this velocity. However, strong spectral 
lines at these wavelengths are so rare that this seems highly un- 
likely. At the moment, we can only speculate that the anomalous 
hyperfine ratio may be due to some non-LTE excitation mecha- 
nism, possibly caused by hyperfine- selective radiative pumping 
by emission lines of the main isotope at shorter wavelengths. 



4.3. Physical conditions 

The integrated intensity in the three hyperfine components adds 
up to 52 K km/s, which - in the high density limit (n > 
3000 cm -3 ) - converts into a colu mn density N( l3 C + ) = 2.6 x 
10 17 cm -2 (Crawfo rd et al.l Il985l) . Assuming standard abun- 
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tiofi 



dance ratiofl this corresponds to yV( 12 C + ) = 1.6 x 10 9 cm 



r- 2 
and a total hydrogen column density of N(H) = 1.6 x 10 23 cm -2 
in the gas emitting the [Cn] lines, if all the carbon is ionized. 

Owing to the foreground absorption, the intrinsic [ 12 Cn] 
emission is not constrained well. The background emission line 
component requires temperatures of several 100 K, with 165 K 




Fig. 2. [ 12 Cn] channel maps of 1 km/s wide velocity channels 
spaced by 1 km/s. Panels are labeled by their center velocity. 
The markers denote the positions of the clumps FIR2 to FIR6 
(lMezgeretal.1119881) . Contour levels are from 10 K km/s to 150 
K km/s in steps of 20 K km/s. 



Assumed abundances: 12 C/ 1J C = 60; 12 C/H = 10 



NGC 2024 
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Fig. 3. Spectrum averaged over the area outlined in Fig. Q~|with 
marks indicating the predicted ve locities of the [ 13 Cn] hyperfine 
components (ICooksy et al. II 19861) . 
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Fig. 4. Spectrum of Fig.[3]replotted on a common vlsr scale and 
overlaid by the model fit (see text and Table [T]). The dashed line 
marks the center velocity of the [ 12 Cn] background component. 
Horizontal bars indicate the uncertainty in the rest frequencies. 
The [ 12 Cn] has been divided by ten. 



being a hard lower limi£|. For a gas temperature of e.g. 400 K, 
the [ 12 Cn] line would be almost optically thick (r « 2.4) with an 
intrinsic line temperature of 7mb([ 12 Cii] background) « 325 K. 
Standard spheric al photon-domina ted region (PDR) models 
(Storzer et al. 1996; Rollig et al. 2006) of even the most extreme 
single clumps fall short by about one order of magnitude of re- 
producing this high line intensity and large column density of 
ionized carbon. It is plausible that the edge-on geometry of the 
ionization front may significantly boost the total column density. 
Thus, an equivalent of ten extreme clumps would be required to 
model the observed emission. Similarly, in a clumpy PDR model 
dCubick et al .1120081) . where the amount of material exposed to 
UV radiation is drastically higher, the predicted [Cn] emission is 
much stronger. For instance, a typical modefl with a total mass 
within the beam of 1 M , a radiation field of x - 10 3 an d an 
average density of 10 6 cm -3 can produce the observed [ 13 Cn] in- 
tensity. These values correspond to a hydrogen column density 
of ^10 23 cm" 2 , in agreement with the above estimate. 



3 Owing to the Rayleigh- Jeans correction, the peak T MB of 125 K 
(Fig.|4]) requires gas temperatures of >165 K 

4 The single line data are insufficient to fully constrain a model 



Table 1. Parameters of the fit to the spectrum of Fig. [4] £vlsr 
denotes the velocity offset of the [ 13 Cn] hyperfine satellites to the 
[ 12 Cn] background c omponent. Bracket symbols: []: fixed value; 
{}: literature value (iCooksy et al.lll986b : (): weakly constrained 
value. 



~FCn] 

;s. absorb. 



1 ->0 



2^ 1 



1 -> 1 



?MB [K] 



[400] 



-352.5 



4.3 



8.2 



(1.2) 



vlsr [km/s] 


10.26 


10.20 


-54.9 


21.9 


(72.7) 


Avlsr [km/s] 


3.8 


2.7 


[3.8] 


[3.8] 


[3.8] 


<5v LS r [km/s] 


— 


— 


-65.2 


11.7 


(62.4) 




— 


— 


{-65.2} 


{11.2} 


{63.2} 


rel. intens. 


— 


— 


0.31 


0.60 


(0.09) 




— 


— 


{0.356} 


{0.444} 


{0.200} 



With the absorption dip reaching down to 53 K (Fig. [4]), the 
temperature in the foreground gas is limited to T?q < 90 K0. 
The optical depth derived for the foreground gas is tfg > 1, 
where the exact value depends on the assumptions made about 
the temperatures. If we accept an intrinsic main beam brightness 
temperature of 325 K for the background component, the col- 
umn density of ionized carbon in the absorbing layer needs to be 
N(C + ) > 10 18 cm" 2 for any gas temperature above 40 K. 



5. Conclusions 

We have observed both [ 12 Cn] and [ 13 Cn] emission toward NGC 
2024. The spatial distribution of the emission follows the 8 fim 
continuum emission, tracing the interface between the Hn region 
and the molecular cloud. The isotopic line emission detected 
from an area near the ionization front implies that: 

- The double-peaked spectrum of the [ 12 Cn] line is due to 
strong background emission absorbed by cooler foreground 
gas at the same velocity. 

- The kinematic signature of the ionized gas differs from that 
of the molecular gas as traced by rare isotopes of CO. 

- The total hydrogen column density of warm gas at sev- 
eral hundred Kelvin traced by the [Cn] emission amounts to 
>10 23 cm" 2 , suggesting that the interface is highly clumpy. 

- The strong foreground absorption implies a hydrogen col- 
umn density >10 22 cm" 2 of gas at temperatures below 100 K. 

- The measured intensity ratio of the hyperfine components 
does not match the theoretical value, suggesting that there is 
some contribution from a non-LTE excitation effect 
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